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We report a high-resolution angle-resolved photo-emission spectroscopy study of the evolution of the elec¬ 
tronic structure of FeSei_ccScc single crystals. Isovalent S substitution onto the Se site constitutes a chemical 
pressure which subtly modifies the electronic structure of FeSe at high temperatures and induces a suppression 
of the tetragonal-symmetry-breaking structural transition temperature from 87 K to 58 K for a; = 0.15. With 
increasing S substitution, we find smaller splitting between bands with dyz and dxz orbital character and weaker 
anisotropic distortions of the low temperature Fermi surfaces. These effects evolve systematically as a function 
of both S substitution and temperature, providing strong evidence that an orbital ordering is the underlying order 
parameter of the structural transition in FeSei-ccS^. Finally, we detect the small inner hole pocket for x=0.12, 
which is pushed below the Fermi level in the orbitally-ordered low temperature Fermi surface of FeSe. 


The nature and origin of the ordered phases found in 
proximity to high temperature superconductivity relate di¬ 
rectly to the mechanism of unconventional superconductiv- 
ity m In the iron-based superconductors, the parent non¬ 
superconducting phases typically exhibit a magnetic stripe or¬ 
dering, which breaks four-fold symmetry. However in cer¬ 
tain regions of the phase diagram, the four-fold symmetry 
of the lattice is broken at a higher temperature than the 
onset of magnetism. One explanation is that this so-called 
nematic phase nevertheless originates from magnetic interac¬ 
tions, which in some models may pick out a particular di¬ 
rection to break fourfold symmetry without long-range mag¬ 
netic ordering (H. Alternatively, it has been suggested that 
Ts corresponds to an orbital instability breaking the degen¬ 
eracy of the Fe dxz and dyz orbitals, distinct from magnetic 
order 1213 . FeSe provides a fascinating test case for these 
ideas, since it undergoes a structural transition at Ts=87 K but 
does not magnetically order (41 and has been the subject of 
intense theoretical study Ed. Detailed studies of the elec¬ 
tronic structure in the tetragonal-symmetry-broken phase la- 
[13 have found a prominent splitting of bands with dxz and 
dyz orbital character, which also manifests in strong distor¬ 
tions of the Fermi surface. However, this orbital ordering has 
been found to be sensitive to both physical pressure dlMl 
and substrate-induced strain in multilayer thin film samples 
Ellis. Chemical pressure, as in isovalently substituted 
FeSei-a^Sa;, provides another tuning parameter which can be 
used to study how the electronic structure evolves as orbital 
order is suppressed, helping to establish the nature of the or¬ 
der parameter associated with Tg , and determining to what ex¬ 
tent the orbital fluctuations found in proximity to the ordered 
phase may play a role in the unconventional superconductiv¬ 
ity. 

In this paper we report a high-resolution angle-resolved 
photo-emission spectroscopy (ARPES) study of the electronic 
structure of S-substituted ^^-FeSe single crystals. In the high- 
temperature tetragonal phase, we find that the sizes of the 
Fermi surfaces and the Fermi velocities both increase with S 
substitution. The suppression of the structural transition Tg 


by S substitution correlates with a reduction of the electronic 
anisotropy and the splitting of bands with dxz and dyz orbital 
characters, which onset at Tg. This provides strong evidence 
that the instability towards orbital ordering drives the fourfold 
symmetry-breaking transition. Finally, we show that the in¬ 
ner hole pocket which is pushed below the Fermi level in the 
orbitally-ordered low temperature Fermi surface of FeSe (H 
does cross the Fermi level in FeSei-^^S^^ for x=0.12. 

Experimental Details. Samples were grown by the 
KCI/AICI3 chemical vapour transport method (TTl [T^ . 
ARPES measurements were performed at the 105 beamline at 
the Diamond Light Source, UK. Single crystal samples with 
sizes ~ 1x1 mm were cleaved in-situ below 15 K at a 

pressure lower than 2x 10“^^ mbar. Samples were carefully 
oriented with the scattering plane along M-F-M (Fe-Fe bond 
direction). Measurements were performed using linearly- 
polarised synchrotron light at 20-80 eV, employing a Scienta 
R4000 hemispherical electron energy analyser. Band structure 
calculations were performed in Wien2k |[T9l using the GGA 
approximation (201. 

High temperature ARPES data. Fig. shows a compari¬ 
son of ARPES data of FeSei-a^S^^ single crystals for x=0 and 
x=0.15 at temperatures above Tg. The quality of the ARPES 
spectra of the substituted crystals remains high, since the Fe 
3d orbitals which dominate the low-energy electronic struc¬ 
ture are not directly affected by the substitutional disorder 
outside of the Fe plane. The comparison between the ARPES 
spectra and the density functional theory calculations (DFT) 
in Fig. [2 reveals several important differences, as has been 
previously reported EHni. Firstly, ah observed pockets are 
substantially smaller than DFT predictions 0. Secondly, the 
predicted 7 hole band with dxy character (Fig.^c)) is found 
to not cross the Fermi level but remains ^50 meV below it 
(Fig. [Jd)) with a particularly large band renormalisation of 
^8 compared to the other bands EKH. Thirdly, the outer 
electron band 6 with dxy character is not observed in experi- 
ments (Fig-I^g)). 

The effect of substituting smaller S ions onto the Se site 
results in a lattice contraction El and intuitively might re- 
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FIG. 1. (color online) High temperature electronic structure of 
FeSei-ccScc. ARPES data for a) x=0 and b) x=0.15 for the hole-like 
bands at the Z point (measured using a photon energy of 23 eV). c) 
Band dispersions from a DFT calculation at the Z point, compared 
to d) the extracted peak positions (from fits to MDCs) for x=0 and 
x=0.15. e-h) ARPES data, DET calculation and extracted band po¬ 
sitions for the electron pockets at the M point (at 37 eV). The pre¬ 
dicted outer electron band S with dxy orbital character is not ob¬ 
served. Bands are labelled a — e, as in previous studies Una. 


suit in a greater orbital overlap, increasing the bandwidth. In 
order to perform a quantitative comparison, we extract the 
band positions from constrained fits to the momentum distri¬ 
bution curves (MDCs), focusing on the high-symmetry cuts 
through the Z and M points, as shown in Fig. [2d,h) for x=0 
and x=0.15 and summarized in Table I in the Supplementary 
Material (SM). We find two systematic features as a function 
of S substitution; the Fermi surfaces generally increase in size, 
and additionally the Fermi velocities (dE/dk at Ep) increase, 
although both only by ^10%. These two trends may explain 
the suppression of T^, since an increase in the size of Fermi 


surfaces is likely to reduce the degree of particle-hole Fermi 
surface nesting (which would be perfect in the limit of point¬ 
like Fermi surfaces), while the increase in Fermi velocities 
corresponds to a reduction in the density of states at the Fermi 
level, reducing the bare susceptibility and thereby suppressing 
the orbital instability. 

Suppression of orbital ordering with S substitution. In 
Fig.[^a-h) we show how the electron pockets observed in a cut 
through the M point evolve with S substitution, above and be¬ 
low Tq. At high temperatures, only the e band (with d^zjdyz 
orbital character and defined in Fig[^ is observed (221. A 
prominent feature of ARPES studies of FeSe in both bulk 
crystals and multilayer thin films is a splitting of intensity at 
the M point below T^, which we denote Am- This splitting, 
reaching Am = 49(2) meV for bulk FeSe, corresponds to a 
splitting of the bands with dyz and dxz orbital characters at the 
M point in the Brillouin zone EHm, which are degenerate in 
the high-temperature fourfold-symmetric phase. The value of 
Am is much larger than what would be expected from DFT 
calculations simply taking into account the small orthorhom¬ 
bic distortion of the lattice M, which was taken to be one 
indication that the structural transition at is primarily an 
electronic instability, and more specifically due to orbital or¬ 
dering. Therefore, the splitting Am, extracted from the posi¬ 
tions of peaks in the energy distribution curve (EDC) at the M 
point (Fig. [^k)), is one experimental proxy order parameter 
for the orbital ordering which can be studied as a function of 
temperature and S substitution. 

The orbital ordering, which can alternatively be described 
as the development of an orbital polarisation An = n^z —^yz, 
also manifests as a distortion of the Fermi surface. As shown 
in the low-temperature Fermi surface map in Fig. |^)), sec¬ 
tions of the electron pockets with dyz orbital character con¬ 
tract while dxz sections expand, leaving an elongated and 
anisotropic Fermi surface, which is observed to form a cross¬ 
shape due to effect of sample twinning. Experimentally, we 
can extract the Fermi /ci?-vector of the dyz portion of the elec¬ 
tron pocket from peak fits of the MDC at Ep (see Fig.|^e-h)) 
and follow the Fermi surface distortion as a function of tem¬ 
perature, as shown in Fig. [^1)). This provides another proxy 
order parameter for the orbital ordering, which we define as 
Akp = kp — kpQ n, where kp{) is the Fermi /c-vector just 
above Tq in the tetragonal phase. 

Figs. [^m,n) show how these two proxy order parameters 
for the orbital ordering, the band splitting Am and Fermi sur¬ 
face distortion Akp, both evolve systematically as a func¬ 
tion of temperature and S substitution. The onset of both 
parameters coincides with (as determined separately from 
anomalies in resistivity measurements (23l) and both follow 
an order-parameter-like temperature dependence, as shown in 
Fig.|^m,n). The fact that both order parameters change sys¬ 
tematically with S substitution is strong evidence that the re¬ 
duction of orbital ordering is directly related to the suppres¬ 
sion of Ts in ¥qSci-x^x (see also Fig.|^. 

The emergence of the inner hole band for x=0.12 at low 
temperatures. We now consider the evolution of the hole 
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FIG. 2. (color online) Suppression of orbital ordering with S substitution on the electron pockets, a-h) ARPES data for a cut through 
the M point (F-M-F direction) at low and high temperatures for FeSei-ccScc with a; = 0 — 0.15 showing the prominent splitting of intensity 
below Tg. i) High and j) low-temperature Fermi surface maps for a;=0.15. Cartoon insets show schematic Fermi surface above Tg (top), 
which is distorted below Tg (middle), and the experimental case of twinned crystals where the dxy electron band is not observed (bottom), k) 
Comparison of low temperature EDCs at the M point showing the reduction of the splitting parameter. Am, with S substitution. 1) Temperature 
dependences of the Fermi kp vector of the dyz portion of the Fermi surface, indicating a distortion A/cf which onsets at Tg. The temperature 
dependence of the extracted order parameters of the orbital ordering, m) Am and n)AA;F, as a function of S substitution, as discussed in the 
main text. Solid lines are guides to the eye. 


pockets with S substitution at low temperatures, as shown in 
Fig. [3] Similar to the electron pockets, the degree of Fermi 
surface anisotropy on the hole pockets reduces with S sub¬ 
stitution, as indicated by the reduced splitting of the two kp 
values of the outer hole band, a, as shown in Fig. [^a-f) and 
Table I in SM. In all samples at high temperature the inner [3 
hole band forms a tiny 3D pocket around Z (see Fig. Ea)), 
where a splitting of ^20 meV (H between the a and [3 hole 
bands occurs due to spin-orbit coupling only 1^ . At low 
temperature an extra splitting associated with the orbital or¬ 
dering 1^ pushes the [3 band below the Fermi surface ISlfT^ 
in FeSe, as can be seen by the lack of features in the Fermi 
level MDC for FeSe in Fig. [^a). However, as the orbital or¬ 
dering is reduced with S substitution, this inner pocket can ex¬ 
ist at low temperatures once more, and is clearly observed at 
X = 0.12 (Fig.j^c)). The photon-energy dependence shown 
in Fig. I^g) (which is equivalent to a kz dispersion) shows 
that it forms a 3D pocket around Z, not extending through¬ 
out the Brillouin zone, and also seen in the detailed map of 
the Fermi surface around the Z point in Fig.j^h). This pocket 
will also break four-fold symmetry, as indicated by the inset to 
Fig.j^h). This new hole pocket could appear as an additional 
low frequency in quantum oscillations measurements ll23]| in 
addition to those corresponding to three different bands in pre¬ 
vious reports 

Fig. Qa) shows a proposed phase diagram of FeSei-a^S^^ 
based on the experimental values of the structural and su¬ 
perconducting transition temperatures extracted from our re¬ 
sistivity measurements (to be presented elsewhere 1 ^ ) and 
Refs. lEH [27l. With S substitution there is a strong sup¬ 
pression of the structural transition, whereas the Tc of the 


superconducting state remains relatively constant for low x. 
As shown in Fig. |2b,c), both the low temperature Am and 
Akp/kpo extracted from our ARPES data reduce with x, cor¬ 
relating with Ts in the phase diagram. A similar trend was 
also found in thin films EEl, where Am and Tg are found 
to increase for thinner samples which are more strained. Thus, 
the systematic change of the degree of orbital ordering with S 
substitution, as measured by our two experimental proxy order 
parameters, indicates that the true electronic order parameter 
associated with the structural transition in FeSei-a^S^^ is the 
orbital polarisation An = n^z—^yz- Therefore, the transition 
at Tg corresponds to an instability of the electronic structure 
of this iron-based superconductor in the charge/orbital chan¬ 
nel, which is distinct from the stripe spin density wave phase 
observed in other systems such as BaFe 2 As 2 . 

Studies of FeSe under applied physical pressure have found 
that the structural transition Tg is initially suppressed in a 
similar manner to the chemical pressure effect of S substitu¬ 
tion |[T3l[T4]|, although the magnetic phase (131 [28l and high 
temperature superconductivity 1291 found under higher phys¬ 
ical pressures have not been reported to occur in FeSei-a^S^^ 
(211123 (at ambient pressure ||30l). NMR studies found that 
spin fiuctuations in FeSe increase in magnitude with physical 
pressure EH while Tg is suppressed (OKTdl, suggesting that 
magnetic interactions are not crucial for the structural insta¬ 
bility. On the other hand, our ARPES results show a direct 
correlation between the decrease in Tg and the suppression 
of orbital ordering under chemical pressure (Fig. |^. How¬ 
ever, the fact that the superconductivity in FeSei-a^S^^ seems 
to be almost independent of the orbital ordering and there is 
no pronounced superconducting dome in the absence of any 
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FIG. 3. (color online) Emergence of the small inner hole band 
at low temperatures for x=0.12. a-c) ARPES data for a high- 
symmetry cut through the centre of the Brillouin zone at 23 eV (near 
the Z point) at T J^ll K for x=0-0.12. Solid black lines are the 
MDCs at the Fermi level, indicating the splitting of the outer hole 
band a and the emergence of the inner /3 band for x = 0.12. d-f) 
As above, but for a cut at the F point (37 eV). g) Photon energy- 
dependence of the Fermi level MDC at x = 0.12, indicating that the 
inner hole band forms a 3D pocket around the Z point, as also shown 
by the Fermi surface maps at h) the Z (23 eV) and i) F point (37 eV). 
Cartoon insets show how the apparent splitting of bands arises from 
Fermi surface distortion and sample twinning. 


magnetic phase (Fig. Qa)) counts against orbital fluctuations 
providing a major component of the unconventional pairing 
mechanism. The small size of all Fermi surfaces and the 
absence of the inner hole pocket in particular gives FeSe a 
very different electronic structure to other Fe-based supercon¬ 
ductors. However the emergence of the inner hole pocket in 
FeSei-a^Sa; could have a signiflcant impact on the suscepti¬ 
bility and pairing interactions, and it is also likely to emerge 
when orbital ordering is suppressed under physical pressure in 
FeSe. The exact nature of the orbital ordering is likely to be 
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FIG. 4. a) (color online)Proposed phase diagram of FeSei_ccScc. 
Ts and Tc values based on our resistivity measurements 1^ . com¬ 
pared to literature values of Tc reported in Ref. ED (resistivity mid¬ 
point) and Ref. |27l (heat capacity). 00 and SC refer to orbital or¬ 
der and the superconducting phase respectively. b,c) Evolution of the 
band splitting. Am, and Fermi surface deformation, Akr/kpo, with 
S substitution at low temperature, demonstrating that the magnitude 
of the orbital ordering is directly correlated with the suppression of 
Ts. 


more complex than simply a uniform site-centred ferro-orbital 
ordering between dxz and dyz orbitals ||5l[T2l[24l[32l, and fu¬ 
ture studies using strain to detwin samples may shed light on 
the details of the orbital mechanism. An alternative viewpoint 
is that the Fermi surface deformations which are suppressed 
with S substitution could be interpreted as a manifestation of 
a d-wave Pomeranchuk instability, which is predicted to be 
weakened as the size of Fermi surface increases (331. 

In summary, our high-resolution ARPES study of 
FeSei_a^Sa^ single crystals has established that the chemical 
pressure increases the size of the Fermi surface pockets and 
Fermi velocities at high temperature. Most importantly, we 
have shown the decrease of Tg is associated with a suppres¬ 
sion of orbital order with S substitution, by observing a re¬ 
duced splitting of dxz and dyz orbitals and weaker Fermi sur¬ 
face anisotropy. These results demonstrate that this orbital or¬ 
dering is the true order parameter of the tetragonal-symmetry¬ 
breaking transition at Tg in FeSe. Finally, these subtle changes 
to the electronic structure lead to the appearance of the inner 
hole band in the low temperature Fermi surface. FeSei-a^S^^ 
is emerging as a key system to help resolve long-standing con¬ 
troversies related to role of spin and orbital degrees of freedom 
in the iron-based superconductors. 
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